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This  study  is  basically  composed  of  two  sections.  In  the  first  section,  a  CFD  analysis  is  used  to  provide 
a  better  insight  to  molten  carbonate  fuel  cell  operation  and  performance  characteristics  at  very  high 
current  densities.  Therefore,  a  mathematical  model  is  developed  by  employing  mass  and  momentum 
conservation,  electrochemical  reaction  mechanisms  and  electric  charges.  The  model  results  are  then 
compared  with  the  available  data  for  an  MCFC  unit,  and  a  good  agreement  is  observed.  In  addition,  the 
model  is  applied  to  predict  the  unit  cell  behaviour  at  various  operating  pressures,  temperatures,  and 
cathode  gas  stoichiometric  ratios.  In  the  second  section,  a  thermodynamic  model  is  utilized  to  examine 
energy  efficiency,  exergy  efficiency  and  entropy  generation  of  the  MCFC.  At  low  current  densities,  no 
considerable  difference  in  output  voltage  and  power  is  observed;  however,  for  greater  values  of  current 
densities,  the  difference  is  not  negligible.  If  the  molten  carbonate  fuel  cell  is  to  operate  at  current  densities 
smaller  than  2500  Am-2,  there  is  no  point  to  pressurize  the  system.  If  the  fuel  cell  operates  at  pressures 
greater  than  atmospheric  pressure,  the  unit  cell  cost  could  be  minimized.  In  addition,  various  partial 
pressure  ratios  at  the  cathode  side  demonstrated  nearly  the  same  effect  on  the  performance  of  the  fuel  cell. 
With  a  60  K  change  in  operating  temperature,  almost  10%  improvement  in  energy  and  exergy  efficiencies 
is  obtained.  Both  efficiencies  initially  increase  at  lower  current  densities  and  then  reach  their  maximum 
values  and  ultimately  decrease  with  the  increase  of  current  density.  By  elevating  the  pressure,  both 
energy  and  exergy  efficiencies  of  the  cell  enhance.  In  addition,  higher  operating  pressure  and  temperature 
decrease  the  unit  cell  entropy  generation. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Molten  carbonate  fuel  cells  (MCFCs)  are  electro-chemical  energy 
conversion  devices  which  exhibit  low  emissions,  high  energy 
conversion  efficiency  and  flexibility  regarding  fuel  type,  repre¬ 
senting  a  viable  option  as  future  low-carbon  energy  devices. 
High-temperature  molten  carbonate  fuel  cells  are  particularly 
suited  for  the  stationary  co-generation  of  electrical  power  and  heat, 
and  distributed  energy  supply,  which  enables  the  development  and 
use  of  this  technology  independently  from  the  establishment  of  a 
hydrogen  infrastructure  [1]. 

Despite  the  fact  that  MCFCs  have  been  intensively  investigated 
over  the  past  few  decades,  the  major  commercialization  challenges 
of  molten  carbonate  fuel  cells  such  as  corrosion,  short  lifetime,  rela¬ 
tively  slow  oxygen  reduction  reaction,  low  power  density,  and  high 
cost,  are  still  concerned  [2].  In  general,  the  goals  of  most  ongoing 


*  Corresponding  author.  Tel.:  +1  647  969  5880;  fax:  +1  905  721  3370. 
E-mail  addresses:  Masoud.Youseframandi@uoit.ca  (M.Y.  Ramandi), 
Ibrahim.Dincer@uoit.ca  (I.  Dincer). 

0378-7753 /$  -  see  front  matter  ©  2011  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 1.06.040 


investigations  are  to  lessen  the  cost,  extend  the  cell  life  with  a  stable 
performance  and  make  the  system  more  efficient. 

The  study  of  performance  characteristics  in  MCFCs  is  one  of  the 
fundamental  subjects  of  electrochemistry  and  has  always  been  one 
of  the  major  interests  for  researchers.  Many  studies  can  be  found 
in  literatures  which  have  experimentally  investigated  performance 
of  molten  carbonate  fuel  cells  (e.g.  [3-6]). 

In  addition,  some  researches  [7-9]  are  performed  based  on 
a  derived  ‘formula  for  MCFC  performances’.  On  the  other  hand, 
there  are  also  some  other  studies  which  employ  PDEs,  developed 
by  using  first  principles.  In  literatures,  polarization  effects,  elec¬ 
tric  potentials  and  over-potentials  are  treated  in  various  ways, 
mostly  using  experimentally  derived  correlations  [10-12].  Previ¬ 
ously,  the  volume-averaging  of  variables  in  the  three  phases  was 
used  as  an  investigational  technique  to  develop  a  few  polariza¬ 
tion  models  [13-15].  Based  on  this  technique,  Subramanian  et  al. 
[16]  developed  a  one  dimensional,  steady-state  model  to  inves¬ 
tigate  the  contribution  of  different  cell  components  to  losses  in 
cell  performance  and  showed  that  at  low  current  densities,  the 
electrolyte  matrix  contributes  to  the  major  fraction  of  potential 
losses.  Xu  et  al.  [17]  used  lattice  Boltzmann  method  and  applied 
Brinkman-Forchheimer-extended  Darcy  equations  together  with  a 
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Nomenclature 

a  constant  coefficient  of  enthalpy  and  entropy  equa¬ 

tions 

Av  reaction  surface  area  density  (m2  m-3 ) 

A  cell  active  area  (m2 ) 

C  molar  concentration  (mol  m-3 ) 

cp  specific  heat  (J  kg-1  K-1 ) 

D  mass  diffusivity  of  species  (m2  s-1 ) 

F  Faraday’s  constant,  96,485  (C  mol-1 ) 

h  enthalpy  (J  mol-1 ) 

i  current  density  (A  nrr 2 ) 

Jo  reference  exchange  current  density  (A  m-3 ) 

J  current  density  (A  m-2 ) 

/  volumetric  current  density  (A  nrr3 ) 

I<  permeability  (m2) 

L  length  scale  (m) 

m  mass  flux  (kg  m-2  s-1 ) 

M  molecular  weight  (g  mol-1) 

n  number  of  electrons 

h  molar  flow  rate  (mol  s-1) 

P  static  pressure  (atm) 

Q  transferred  heat  rate  (W) 

R  universal  gas  constant  8.314  (J  mol-1  K-1 ) 

s  entropy  (J  mol-1  K-1 ) 

S  source  terms 

t  time  (s) 

T  temperature  (I<) 

u  gas  velocity  (ms-1) 

W  width  scale  (m) 

W  Power  output  (W) 

x  x-coordinate 

X  molar  fraction 

Greek  letters 

a  transfer  coefficient 

P  reaction  order 

s  porosity 

ri  over-potential  (V),  energy  efficiency  (%) 

fi  dynamic  viscosity  (kg  m-1  s-1 ) 

v  species  stoichiometric  coefficient  of  the  reaction 

0  electric  potential  (V) 

exergy  efficiency  (%) 

E  exergy  (J  mol-1 ) 

p  density  (kg  m-3) 

a  electric  conductivity  (S  m-1 ) 

Subscripts  and  superscripts 
a  anode 

age  anode  gas  channel 

c  cathode/carbonate  ion 

ege  cathode  gas  channel 

chem  chemical 

d  destroyed 

e  electronic 

g  gas  phase 

gen  generation 

in  inlet 

i  ith  component 

j  jth  species 

m  mass  equation 

out  outlet 

r  reversible,  heat  loss  coefficient 

s  solid 


u  momentum  equation 

0  standard  conditions 

0e  electronic  charge  equation 

0c  carbonate  ion  charge  equation 

eff  effective 

ref  reference  state 


reaction-diffusion  equation  with  several  reasonable  assumptions 
to  generate  a  more  realistic  voltage-current  density  curve.  Kim 
et  al.  [18]  also  studied  the  effect  of  water-gas  shift  reaction  on  cell 
performance.  They  figured  out  that  even  though  the  cell  voltage 
and  power  are  slightly  higher  when  the  water-gas  shift  reaction  is 
included,  however  it  is  very  small.  In  stack  level  models,  Yoshiba 
et  al.  [19]  incorporated  numerical  simulation  to  compare  the  stack 
performance  of  various  gas-flow  types.  Furthermore,  an  interest¬ 
ing  study  was  carried  out  by  [20]  to  determine  the  behaviour  of  a 
unit  MCFC  at  the  beginning  of  the  operation.  It  is  worth  to  know 
that  the  time  to  achieve  a  steady-state  was  about  half  a  second.  In  a 
study  [21  ]  for  a  multi-channel  unit  cell,  it  was  found  that  the  non¬ 
uniformity  in  fuel  flow  through  channels  deteriorates  the  fuel  cell 
performance  dramatically. 

On  the  other  hand,  thermodynamic  analysis  has  been  used  as 
a  viable  tool  to  investigate  the  energy  and  exergy  aspects  of  the 
cogeneration  systems.  In  fact,  an  MCFC  permits  the  recovery  of 
waste  heat,  which  can  be  used  in  the  production  of  steam,  hot 
or  cold  water,  or  hot  or  cold  air,  depending  on  the  associated 
recuperation  equipment  [22].  For  instance,  employing  a  bottom¬ 
ing  cycle  consisting  of  a  turbine  will  help  to  supply  the  necessary 
power  required  for  the  compressors.  [23].  As  an  effort,  Silveira  et 
al.  [22],  implemented  the  energy,  exergy  analysis  for  a  cogenera¬ 
tion  system  and  accomplished  a  global  efficiency  or  fuel  utilization 
efficiency  of  86%.  Their  analysis  showed  that  the  exergy  loss  in  the 
fuel  cell  unit  is  significant.  Varbanov  et  al.  [24]  offered  a  broader 
view  of  the  concept  of  power  generation  combined  cycle  by  com¬ 
bining  an  MCFC,  rather  than  a  gas  turbine,  with  a  steam  turbine. 
Their  results  showed  that  the  inherently  high  power  generation 
efficiency  of  the  MCFC  (46.38%  in  this  case)  can  be  significantly 
increased,  up  to  nearly  70%.  Kang  et  al.  [25]  simulated  an  exter¬ 
nally  reformed  MCFC  system  and  analyzed  the  effects  of  system 
configuration  and  operating  conditions  on  the  system  efficiency 
for  a  lOOkW  MCFC  system.  Rashidi  et  al.  [26]  performed  energy 
and  exergy  analyses  of  a  molten  carbonate  fuel  cell  hybrid  sys¬ 
tem.  An  overall  energy  efficiency  of  57.4%,  exergy  efficiency  of 
56.2%,  bottoming  cycle  energy  efficiency  of  24.7%  and  stack  energy 
efficiency  of  43.4%  were  achieved,  accordingly.  Musa  et  al.  [27]  pre¬ 
sented  thermodynamic  models  for  the  internally  and  externally 
reformed  MCFC  and  employed  them  in  an  ASPEN  Customer  Mod¬ 
eller.  The  results  indicated  that  the  operating  temperature  has  more 
effect  on  the  cell  voltage  of  IR-MCFC  system  compared  to  ER-MCFC 
system. 

Considering  all  previous  studies,  it  is  noticeable  that  the  high 
cost  of  MCFCs  is  still  an  open  challenge  for  commercialization. 
Operation  at  high  current  densities  would  be  an  option  to  lessen 
cell  cost.  Most  of  the  previous  studies  have  focused  on  the  reg¬ 
ular  operating  current  densities  (e.g.  0.1-0.25  A  cm-2).  The  scope 
of  this  research  is  to  provide  a  better  insight  to  molten  car¬ 
bonate  fuel  cell  operation  and  performance  characteristics  at 
very  high  current  densities  (up  to  0.6  A  cm-2)  in  order  to  reduce 
the  molten  carbonate  fuel  cells  size,  weight  and  consequently 
cost,  it  is  essential  to  develop  fuel  cells  with  higher  power  den¬ 
sities.  This  is  achieved  by  strictly  applying  balances  of  mass, 
momentum,  species  and  charges  and  avoiding  highly  algebraic 
empirical  correlations.  Literatures  usually  present  a  polarization 
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Fig.  1.  The  physical  domain  of  the  simulated  MCFC. 


curve  which  is  derived  by  applying  experimental  correlations  for 
electric  potentials  and  over-potentials.  However,  in  this  study,  a 
more  complicated  mathematical  model,  involving  several  coupled 
partial  differential  equations  (PDEs)  is  used.  Each  phenomenon, 
even  electric  charge  transfers  have  a  separate  PDE  which  needs 
to  be  solved  using  computational  fluid  dynamics  (CFD).  With 
the  intention  of  prevent  using  highly  empirical  correlations  for 
electrodes  and  electrolyte  resistance,  a  diffusion  model  involving 
Butler-Volmer  equation  is  applied  to  solve  the  electronic  and  ionic 
potential  fields.  From  the  thermodynamic  point  of  view,  literatures 
typically  offer  analysis  based  on  the  MCFC  integrated  systems  while 
the  exergy  efficiency  of  the  unit  cell  itself  has  been  reported  rarely. 
Hence,  energy  and  exergy  analyses  are  also  implemented  to  exam¬ 
ine  efficiencies  and  entropy  generation  of  the  unit  cell.  After  all, 
it  is  worthwhile  to  point  out  that  in  the  first  part  of  this  study, 
a  CFD  analysis  is  carried  out  to  obtain  the  characteristics  of  an 
MCFC  and  thermodynamic  analysis  is  put  into  practice  to  study 
the  entropy  generation,  energy  and  exergy  efficiencies  of  the  unit 
cell. 

2.  Problem  formulation 

To  mathematically  describe  a  MCFC,  a  comprehensive  knowl¬ 
edge  of  physical  and  electrochemical  mechanisms  is  required.  The 
proposed  model  for  the  molten  carbonate  fuel  cell,  involves  the 
conservation  of  mass,  momentum,  chemical  species  and  electric 
charges. 

The  physical  domain  of  the  simulated  MCFC  is  shown  in  Fig.  1. 
As  it  may  be  observed,  it  can  be  physically  broken  down  into  five 
individual  zones:  anode  gas  channel  (AGC),  anode,  electrolyte,  cath¬ 
ode  and  cathode  gas  channel  (CGC).  Making  a  theoretically  rigorous 
fuel  cell  model  which  reflects  the  micro/macro-scale  transport  pro¬ 
cesses  is  extremely  challenging  because  of  a  lack  of  experimentally 
evaluated  physical  parameters.  Therefore,  every  numerical  simu¬ 
lation  is  conceived  and  developed  based  on  a  set  of  assumptions. 
The  following  assumptions  are  made  for  the  present  model: 

•  The  gas  flow  in  the  anode  and  cathode  gas  channels  is  laminar 
flow. 

•  The  chemical  species  obey  the  ideal  gas  law  and  are  ideally  mixed. 

•  Water  exists  only  in  gaseous  form. 

•  The  porous  anode  and  cathode  are  homogeneous. 

•  The  effects  of  gravity  are  negligible. 

•  All  reactions  take  place  in  the  triple  boundary  phase. 

•  The  Butler-Volmer  equation  is  employed  for  the  reaction  in  the 
electrodes. 

•  Any  change  in  concentration  of  carbonate  ion,  C032-,  is  neglected. 

•  The  gas  mixture  and  solid  components  of  fuel  cell  are  in  a  thermal 
equilibrium  state. 

•  The  cell  is  operated  with  pure  hydrogen,  hence  no  contamination 
effects  are  considered.  It  is  worth  to  mention  that,  in  this  study, 
fuel  is  considered  to  be  supplied  by  an  external  reformer.  In  this 
reformer,  natural  gas  which  is  previously  mixed  with  steam  and 
preheated  in  a  regenerative  heat  exchanger,  converts  to  hydro¬ 
gen  and  carbon  dioxide.  The  reformed  gas  passes  through  the 


Table  1 

Source  terms  in  conservation  equations. 


Zone 

Sm 

Su 

Si 
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S(j)c 
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0 
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0 
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0 

Anode 

Sh2  +Sco2  +Sh2o 

K~  US 

Sh2  =  - 

J'a 

nF 

J'a 

-J'a 

8 

II 

J'a 

nF 

Electrolyte 

0 

0 

O  gl 

*0 

II 

Ja 

nF 

0 

0 

Cathode 

S  o2  +Sco2 

-e/z  11 

K  US 

So2  =  - 

1  J'c 

~  2  nF 

-J'c 

J'c 

CGC 

0 

0 

II 

J?o 

J'c 

nF 

0 

0 

exchanger  over  and  over  again  to  reduce  the  amount  of  contam¬ 
inants.  Next,  this  mixture  which  has  a  trace  amount  of  methane 
and  carbon  monoxide  is  fed  to  the  anode  side  of  the  MCFC.  It  is 
assumed  that  repeating  the  reforming  process  and  increasing  the 
residual  time  in  the  reformer  have  caused  a  significant  reduction 
in  the  fraction  of  methane  and  carbon  monoxide  and  hence  the 
effect  of  contamination  is  neglected.  This  assumption  was  imple¬ 
mented  by  some  other  researchers  who  have  used  pure  hydrogen 
in  their  simulation  (e.g.  [11,17,21]).  However,  the  contaminant 
effects  will  be  considered  in  our  next  study  which  will  discuss 
both  externally  and  internally  reformed  fuels. 

2.1.  Governing  equations 

In  general,  in  MCFCs  fuel  gas  which  is  a  gaseous  mixture  of 
hydrogen,  water  vapour  and  carbon-dioxide  enters  the  AGC  and 
diffuses  through  the  porous  anode  where  hydrogen  molecules  get 
involved  into  the  hydrogen  oxidation  reaction  (HOR).  During  this 
electrochemical  reaction,  hydrogen  combines  with  carbonate  ions, 
ends  up  to  the  water  vapour  and  carbon  dioxide  generation.  More¬ 
over,  the  released  electrons  migrate  through  an  external  circuit, 
create  electricity  and  return  in  to  the  cell  through  cathode.  On 
the  cathode  side,  a  mixture  of  oxygen,  carbon-dioxide  and  nitro¬ 
gen  enters  the  cathode  gas  channel  (CGC)  and  diffuses  through 
the  porous  cathode,  where  oxygen  reduction  reaction  (ORR)  takes 
place.  Oxygen  is  reduced  to  carbonate  ion  by  combining  with  car¬ 
bon  dioxide  and  the  electrons  coming  from  the  external  circuit.  The 
carbonate  ions  formed  at  the  cathode  move  through  the  electrolyte 
toward  the  anode,  carrying  the  electric  current  and  completing 
the  carbon  dioxide  circuit.  The  underlying  transport  phenomena 
within  the  cell  are  modeled  with  conservation  equations  for  mass, 
momentum,  species  and  electric  charges.  This  section  exhibits  the 
general  form  of  these  equations  applied  to  all  MCFC  zones  and 
the  sink/source  terms  for  each  individual  zone  are  summarized  in 
Table  1. 

The  equation  for  conservation  of  mass,  or  continuity  equation, 
can  be  written  as  follows  [28]: 

9 

+  V  •  (PgUg)  —  Sm  (1) 

where  pg  is  the  density  and  calculated  by 

N 

Pg  =  ^T/W.C,  (2) 

1=1 

In  addition,  Sm  represents  the  volumetric  rate  of  mass  produc¬ 
tion/consumption  which  is  related  to  the  Butler-Volmer  equation 
as  described  at  the  end  of  this  section. 

The  equation  for  conservation  of  momentum  is  [28]: 

^  (“Pg^g)  +  ^  '  (^2  PgUgUg'j  =  Pg  +  ^  •  (/rVUg)  +  SU  (3) 
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where  Su  is  the  source/sink  term  of  the  momentum  equation  which 
is  applicable  just  for  the  porous  electrodes. 

In  fact,  in  this  approach,  the  Darcy  equation  is  utilized  to  model 
porous  media.  This  momentum  sink,  contributes  to  the  pressure 
gradient  in  the  porous  electrodes,  creating  a  pressure  drop  that 
is  proportional  to  the  fluid  velocity.  For  a  homogeneous  porous 
media; 

Su  =  --j^Ug  (4) 

where  I<  is  the  permeability  of  the  porous  media  and  is  determined 
by 


Dj  £ 


150(1_£)3 


(5) 


To  describe  the  chemical  species  transport,  Eq.  (6)  is  considered 
which  is  the  general  form  of  the  conservation  equation  including 
both  convection  and  diffusion  terms  [28]. 


l(sC,)  +  V  ■  (-DfmV C,)  +  V  •  (tigC,)  =  S,  (6) 

The  effective  mass  diffusion  coefficient  for  species  i  in  the  mix¬ 
ture,  Dfjn ,  is  calculated  by 


Deff  = 


1  -X, 


(7) 


where  Djjf  is  the  effective  binary  mass  diffusion  coefficient  of  com¬ 
ponent  i  in  component  j  and  is  calculated  using 


De//  =  Di 


U  T 

1  r 


ref  c1.5 


ref 


(8) 


Before  moving  on  to  the  electric  transport  equations,  it  is  worth 
to  mention  that  in  a  porous  electrode  the  local  electrochemical 
reaction  rate  depends  on  the  driving  force,  the  local  over-potential, 
i.e.,  the  divergence  from  the  equilibrium  potential.  Therefore,  it  is 
preferable  to  use  the  potential  instead  of  the  current  in  the  trans¬ 
port  equation  [29].  Furthermore,  the  current  density  is  a  vector 
quantity,  and  only  the  transverse  component  (normal  to  the  elec¬ 
trode  surface)  is  functional  and  contributes  to  the  power  output  of 
the  cell;  the  lateral  component  only  decreases  the  cell  output  [30]. 

The  charge  conservation  equations  describe  the  electric  current 
in  electrically  conductive  components  and  ionic  current  in  ionic 
conductive  components.  Any  change  in  concentration  of  carbonate 
ion,  C032-,  is  neglected  which  consequently  means  that  the  effect 
of  migration  can  be  neglected  as  well.  Therefore,  Ohms’  law  is  valid 
in  MCFC  zones  for  both  electronic  charge  and  ionic  charge.  Flence, 
the  conservation  of  electronic  charge  can  be  written  as 


0  =  V-(— S1'5fTeV(/>e)-S^  (9) 

and  the  conservation  equations  of  ionic  charge  is 

0  =  V  •  (— e1'5crcV(/>c)  —  S<pc  (10) 

It  should  be  pointed  out  that,  with  the  intention  of  simplifying, 
the  unit  cell  is  considered  to  operate  in  isothermal  condition.  It 
is  worthwhile  to  mention  that  each  transport  phenomena  is  rep¬ 
resented  by  a  general  form  of  conservation  equation  which  takes 
the  porosity  into  consideration.  Clearly,  the  porosity  ranges  from 
zero  to  one  for  various  zones.  Finally,  the  Buttler-Volmer  equation 
which  describes  the  net  withdrawn  current  density  in  fuel  cells  is 
used  to  determine  various  source  terms  of  the  MCFC. 


Ja=Av,aJofa 


exp 


/  Oia,anFl1a 
V  RT 


) 


-  exp 


/  ac,anFria\ 
V  RT  J 


(13) 


J'c=A,,cjfc 


exp 


faa,cnFr]c 
V  RT 


) 


-  exp 


/  ac,cnFrjc\ 
V  RT  J 


(14) 


2.2.  Boundary  conditions 

With  the  purpose  of  completing  the  MCFC  model  formulation, 
specifying  various  boundary  conditions  at  different  positions  are 
critical.  The  boundary  conditions,  for  a  computational  domain  with 
a  single  pair  of  gas  flow  channels,  are  illustrated  in  Fig.  2.  At  the 
anode  gas  channel  inlet  ( Iagc )  and  cathode  gas  channel  inlet  ( Icgc ), 
the  total  mass  flux  and  gas  species  composition  of  the  entering  gas 
flow  are  specified.  Moreover,  the  fluxes  of  electric  and  ionic  charge 
are  considered  to  be  zero.  Additionally,  considering  the  very  large 
aspect  ratio  (length  to  height  ratio)  of  the  gas  channels,  the  flow 
is  assumed  to  be  fully  developed  at  the  anode  gas  channel  outlet 
( Oagc )  and  cathode  gas  channel  outlet  ( Ocgc ).  This  means,  none  of 
the  variable  fluxes  vary  in  the  normal  direction.  Likewise,  the  gas 
pressure  is  specified.  A  no  slip  boundary  condition  is  applied  to  the 
anode  and  cathode  gas  channel  walls  including,  Wjgc,  WRgc. 

The  fluxes  of  electric  and  ionic  charge  are  set  to  zero.  For  all 
other  walls  of  the  gas  channels,  w£agc,  WRagc,  w£cgc,  and  WRcgc, 
the  electronic  potential  is  specified  while,  the  flux  of  the  remaining 
variables  is  set  as  zero.  For  the  no-flux  boundaries,  ZL  and  ZR ,  the 
zero  flux  condition  is  assumed  for  all  variables. 


3.  Numerical  study 

The  Finite  Volume  Method  is  employed  to  simulate  the  opera¬ 
tion  of  the  molten  carbonate  fuel  cells.  The  computational  domain 
is  defined  and  divided  into  a  number  of  control  volumes  using 
ANSYS ICEM  CFD  12.0.1.  To  handle  the  divergence  difficulties,  a  set 
of  under-relaxation  techniques  is  developed.  The  governing  equa¬ 
tions  are  discretized  using  commercial  software,  ANSYS  FLUENT 
12.0.1.  The  C  programming  language  is  used  to  develop  a  code  in 
order  to  add  several  capabilities  to  ANSYS  FLUENT  12.0.1.  The  SIM¬ 
PLE  algorithm  is  selected  for  the  coupling  between  the  pressure 
and  velocity  field.  An  algebraic  multi-grid  (AMG)  method  with  a 
Gauss-Seidel  type  smoother  is  used  to  accelerate  the  convergence. 
A  strict  convergence  criterion  with  a  residual  of  10-8  is  used  for  all 
variables. 


4.  Thermodynamic  analysis 

Thermodynamics  plays  a  critical  role  in  the  analysis  of  processes, 
systems  and  devices  in  which  energy  transfers  and  energy  trans¬ 
formations  take  place  [31].  In  this  regard,  energy  analysis  is  the 
traditional  method  of  assessing  the  way  energy  is  used.  However, 
an  energy  balance  provides  no  information  on  the  degradation  of 
energy  or  resources  and  does  not  quantify  the  usefulness  or  qual¬ 
ity  of  the  various  energy  and  material  streams  flowing  through  a 
system  and  exiting  as  products  and  wastes  [31  ].  On  the  other  hand, 
exergy  analysis  is  a  useful  tool  for  furthering  the  goal  of  more  effi¬ 
cient  energy  use,  as  it  enables  the  determination  of  the  location, 
type  and  true  magnitude  of  energy  wastes  and  losses  in  a  system 
[32].  Exergy  is  defined  as  the  maximum  amount  of  work  which  can 
be  obtained  from  a  system  or  a  flow  of  matter  when  it  is  brought 
reversibly  to  equilibrium  with  the  reference  environment.  Exergy 
analysis  is  based  on  second  law  of  thermodynamics  and  the  con¬ 
cept  of  irreversible  entropy  production.  The  exergy  consumption 
during  a  process  is  proportional  to  the  entropy  production  due  to 
irreversibilities. 
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Fig.  2.  Boundary  conditions  for  the  MCFC  model:  (a) 
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In  order  to  investigate  the  performance  of  the  MCFC  in  thermo¬ 
dynamics  point  of  view,  the  molten  carbonate  fuel  cell  is  treated  as 
a  black  box  with  two  inputs  and  two  output  streams  (Fig.  3)  and 
the  energy  and  exergy  balances  are  applied  for  the  unit  cell.  In  this 
fashion,  energy  and  exergy  efficiencies  of  the  unit  cell  are  defined 
accordingly. 

It  is  worthwhile  to  mention  that  this  research  is  to  present  per¬ 
formance  characteristics  of  a  ‘unit  cell’.  Generally,  the  inlet  streams 
are  taken  from  external  devices  as  they  are.  Therefore,  no  calcula¬ 
tion  is  performed  out  of  the  unit  cell  boundaries,  as  estimating  the 
parasitic  losses  is  not  in  the  scope  of  this  study. 

First  of  all,  for  a  specified  operating  current  density,  components 
molar  usage  and  production  can  be  calculated.  Therefore,  a  molar 
balance  results 
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The  energy  balance  of  the  unit  cell  can  be  written  as  follows: 

N  N  N  N 

^  ^a,inj^a,in  ^  ^cjn^cfn  _  ^  ^a,out j^a,out  ^  ^c,out j^c,out 
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(20) 


Accordingly,  the  exergy  balance  reads 
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(21) 


where  Qmcfc.  Wmcfc  and  Sd  denote  the  heat  transferred  by  MCFC 
walls,  the  output  electrical  power  of  the  MCFC  and  the  destroyed 
exergy  ( T.sgen ),  respectively. 

The  molar  exergy  term  is  evaluated  by  a  summation  over  phys¬ 
ical,  thermal  and  chemical  exergetic  terms.  Therefore 


Si  =  ([hf(T)  -  hf(T0)]  -  T0[s,(T)  -  s,(To)])  +  RT0  ln(X,)  +  S[hem 

(22) 

where  Xz  is  the  molar  fraction  of  component  i  and  S^hem  is  the 
chemical  exergy  which  are  found  from  [33]. 

With  the  purpose  of  evaluating  molar  enthalpy  and  entropy  of 
each  species,  and  assuming  that  the  gases  in  anode  gas  channel  and 
cathode  gas  channel  obey  the  ideal  gas  behaviour,  the  following 
polynomial  equations  fitted  to  the  data,  taken  from  JANAF  table 
used  in  [34]. 

sH“l  +  yT+ +  +  +  1231 

^  =  (fll  In (T)  +  a2T  +  -fT2  +  ^T3  +  ^T4  +  a7)  (24) 

Furthermore,  the  output  power  of  the  unit  cell  is  calculated  as 
follows: 


WmCFC  =  *MCFC  •  ^MCFC  '  ^MCFC  (25) 

where  AMcfc  is  the  active  surface  area  of  MCFC  and  i’mcfc  is  the 
operating  current  density.  The  cell  voltage  is  determined  by 


^MCFC  =  Er  -  1’mCFC (rjact  +  Vconc  +  *]ohm) 


(26) 


Fig.  3.  Schematic  of  the  fuel  cell  and  system  boundaries  studied  in  energy  and  exergy  where  Er  is  the  reversible  Open  circuit  voltage.  Likewise,  Tjact »  *7  cone 

analyses.  and  r]ohm  are  the  activation,  concentration  and  ohmic  impedances, 
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Table  2 

The  structural  parameters  of  the  simulated  MCFC. 


Parameter 

Value 

Anode  gas  channel  height  (mm) 

0.2 

Anode  height  (mm) 

0.07 

Electrolyte  height  (mm) 

1.0 

Cathode  height  (mm) 

0.06 

Cathode  gas  channel  height  (mm) 

0.2 

Cell  length  (mm) 

100 

Cell  width  (mm) 

100 

Porosity  of  anode,  sa 

0.52 

Porosity  of  cathode,  sc 

0.62 

respectively.  These  parameters  are  estimated  by  the  available 
empirical  correlations  in  the  literature  [35]. 
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rjohm  =  0.5  x  10  4  x  exp 

In  addition,  the  heat  transferred  by  MCFC  walls  is 
Qmcfc  =  T(  As  -  Sgen) 
where  Sgen  is  calculated  by 
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Finally,  the  energy  and  exergy  efficiencies  of  the  molten  carbon¬ 
ate  fuel  cell  are  evaluated  by 
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Fig.  4.  Comparison  of  model  prediction  with  the  literature  data  [17]. 
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The  fuel  gas  component  mole  fractions  for  H2,  FI20  and  C02  were 
0.7, 0.2  and  0.1,  respectively.  Furthermore  oxidant  was  assumed  to 
be  consisting  of  02,  C02  and  N2  with  the  mole  fractions  of  0.1 5, 0.30 
and  0.55,  correspondingly.  From  the  results  shown  in  Fig.  5,  it  may 
be  observed  that  for  low  current  densities  (up  to  2000  A  m-2 ),  there 
is  not  a  significant  difference  in  cell  output  voltage  and  power,  but 
for  higher  values  of  current  densities,  the  difference  is  not  negli¬ 
gible.  More  specifically,  a  20  K  temperature  reduction  lowers  the 
output  power  almost  20%.  In  addition,  it  is  clear  that  an  increase  in 
operating  temperature  results  a  higher  output  voltage.  As  it  is  men¬ 
tioned  earlier,  fuel  cells  with  high  current  densities  are  demanded 
to  reduce  the  cost.  But  Fig.  5  shows  that  the  cell  output  voltage 
decreases  as  the  current  density  increases.  Therefore,  an  investiga¬ 
tion  is  required  to  figure  out  what  the  optimum  operating  current 
density  is.  This  could  be  achieved  by  plotting  power  versus  cur¬ 
rent  density  at  different  operating  temperatures  as  shown  in  Fig.  6. 
According  to  this  figure,  the  output  power  grows  with  an  increase 
in  current  density  up  to  the  point  where  the  overall  trend  begins 
to  be  diminished.  This  is  the  point  in  which  power  has  its  highest 
value  and  current  density  also  is  relatively  high.  This  pair  of  current 
density  and  power  could  be  used  for  unit  cell  design  considera- 


5.1.  CFD  results 

The  system  of  conservation  equations  were  solved  after  fix¬ 
ing  some  field  variables  for  various  operating  conditions  using  the 
specifications  described  in  Table  2.  In  order  to  investigate  the  effects 
of  temperature,  pressure,  current  density  and  cathodic  partial  pres¬ 
sure  ratio,  the  voltage-current  density,  power-current  density  and 
overpotential-temperature  curves  were  plotted  at  different  condi¬ 
tions.  In  a  simple  case,  model  prediction  for  voltage-current  density 
curve  was  compared  with  the  literature  data  [17].  Fig.  4  shows 
the  validated  results.  As  may  be  observed  from  this  figure,  a  good 
agreement  was  obtained. 

Next,  for  each  of  the  cases  that  are  studied  in  this  research,  the 
unit  cell  was  assumed  to  be  isothermal  and  operate  at  a  tempera¬ 
ture  of  883,  903  and  923  K.  The  operating  pressure  of  the  cell  was 
set  to  be  1,  2  and  3  atm  and  the  cathodic  stoichiometric  ratios  of 
0.5, 1.0  and  2.0  were  selected. 

Figs.  5  and  6  demonstrate  the  effects  of  temperature  on  polar¬ 
ization  curves  and  power-current  density  curves,  respectively,  for 
three  different  operating  temperatures  (883,  903  and  923  K).  For 
this  case,  the  operating  pressure  was  assumed  to  be  atmospheric. 


Current  Density  (A  m*2) 


Fig.  5.  Effect  of  temperature  on  polarization  curve  (operating  condition:  atmo¬ 
spheric  pressure,  fuel  gas:  0.7  H2/0.2  H2O/0.1  C02  and  0.15  02/0.3  C02/0.5  N2  as 
oxidant). 
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Fig.  6.  Effect  of  temperature  on  output  power  (operating  condition:  atmospheric 
pressure,  fuel  gas:  0.7  H2/0.2  H2O/0.1  C02  and  0.15  02/0.3  C02/0.5  N2  as  oxidant). 

tions.  For  an  operating  temperature  of  883  K,  the  maximum  power 
occurs  at  the  current  density  of  3000  Am-2  while  a  40 K  tempera¬ 
ture  increase  enables  the  system  to  enhance  its  maximum  power 
by  1000  Am-2.  In  fact,  by  raising  the  unit  cell  operating  tempera¬ 
ture,  the  maximum  power  shifts  to  the  right  side  of  the  horizontal 
axis,  which  consequently  results  in  a  higher  current  density. 

The  variations  in  voltage  losses  with  changes  in  temperature 
are  presented  in  Fig.  7.  For  this  case,  three  different  current  den¬ 
sities  (1500,  2000  and  2500  Am-2)  were  chosen.  The  fuel  gas  and 
oxidant  composition  and  cell  operating  pressure  were  assumed  to 
be  the  same  as  the  previous  case.  This  figure  illustrates  the  rea¬ 
son  behind  the  decreasing  tendency  in  cell  output  voltage  as  the 
current  density  becomes  more  intense.  These  losses  are  caused 
by  the  polarization  effects.  The  study  of  polarization  losses  is  one 
of  the  main  subjects  of  electrochemistry  application  in  fuel  cells. 
The  irreversibilities  occurring  in  a  fuel  cell  electrochemical  reac¬ 
tion  degrade  the  unit  cell  output  voltage  in  form  of  over-potentials. 
According  to  Fig.  7,  the  fuel  cell  voltage  losses  drop  as  the  tem¬ 
perature  increases  in  an  almost  linear  fashion  for  all  operating 
current  densities.  In  addition,  at  any  temperature,  higher  polar¬ 
izations  are  caused  by  greater  operating  current  densities.  In  fact, 


Temperature  (K) 

Fig.  7.  Effect  of  temperature  and  current  density  on  total  polarization  (operating 
condition:  atmospheric  pressure,  fuel  gas:  0.7  H2/0.2  H2O/0.1  C02  and  0.15  02/0.3 
C02/0.5  N2  as  oxidant). 


Current  Density  (A  m'2) 

Fig.  8.  Effect  of  pressure  and  current  density  on  cell  power  (operating  condition: 
923  K,  fuel  gas:  0.7  H2/0.2  H2O/0.1  C02  and  0.15  02/0.3  C02/0.5  N2  as  oxidant). 

as  current  density  increases,  the  electrochemical  reaction  rates 
become  enhanced  which,  consequently,  heightens  the  consump¬ 
tion/production  rate  of  chemical  gaseous  species,  surges  the  system 
irreversibilities  and  hence  creates  voltage  losses  when  the  drawn 
current  density  is  increased. 

Fig.  8  presents  the  output  power  versus  current  density  for 
various  operating  pressure  (1  atm,  2  atm  and  3  atm).  Fuel  gas  and 
oxidant  compositions  were  set  to  be  the  same  as  the  previous  case 
studies  and  the  operating  temperature  was  fixed  at  923  K.  As  shown 
in  this  figure,  it  is  obvious  that  for  low  current  densities  (up  to 
2500  Am-2),  there  is  not  a  significant  difference  in  output  power. 
Therefore,  if  the  molten  carbonate  fuel  cell  is  to  operate  at  current 
densities  smaller  than  2500  Am-2,  there  is  no  point  to  pressurize 
the  system  and  accordingly  augment  the  cost.  However,  it  may  be 
observed  from  this  figure  that  for  any  current  density,  the  output 
power  is  enhanced  as  the  operating  pressure  is  increased.  In  addi¬ 
tion,  the  maximum  power  occurs  at  higher  current  densities  for 
greater  pressures.  It  may  be  perceived  that  if  the  fuel  cell  operates 
at  pressures  greater  than  atmospheric  pressure,  the  system  will  be 
capable  of  generating  greater  values  of  power  which,  economically, 
could  lessen  the  costs.  However,  it  should  be  pointed  out  that  pres¬ 
surizing  the  unit  cell  results  in  higher  partial  pressure  of  the  fuel 
and  oxidant  components  which  contribute  in  the  electrochemical 
reactions  and  subsequently,  decrease  the  irreversibilities  leading 
to  a  higher  output  voltage. 

In  practice,  the  benefits  of  pressurized  operation  are  momen¬ 
tous  just  up  to  about  3-4  bar.  However,  existing  manufacturers  do 
not  operate  their  MCFC  systems  at  very  higher  pressures.  Without 
a  doubt  there  are  disadvantages  associated  with  the  higher  operat¬ 
ing  pressures  brought  by  the  system  design  constraints.  In  general, 
it  is  thought  that  pressurizing  the  MCFC  systems  is  not  economical 
for  the  output  power  less  than  1  MW.  There  are  important  power 
costs  involved  in  compressing  the  reactant  gases.  The  costs  involved 
with  operating  at  high  pressure  are  those  connected  with  the  extra 
mass  involved  with  designing  a  high-pressure  system.  Contrast¬ 
ing  the  benefits  of  increased  pressure  are  the  effects  on  unwanted 
side  reactions  such  as  carbon  deposition  (Boudouard  reaction).  The 
problem  of  leaks  from  high-pressure  systems  is  obviously  another 
crucial  concern.  The  MCFC  system  needs  to  be  designed  with  the 
least  possible  chance  of  any  leaks  and  should  be  monitored  for  such 
leaks  on  a  regular  basis.  Apart  from  the  waste  of  gas,  there  is  also  the 
likelihood  of  the  build-up  of  explosive  mixtures  of  hydrogen  and 
oxygen.  To  sum  up,  it  is  crucial  to  select  the  operating  pressure, 
accurately. 
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Fig.  9.  Effect  of  current  density  and  Pc o2  /po2  ratio  on  cell  power  (operating  condi¬ 
tion:  923  K,  fuel  gas:  0.7  h2/0.2  h2o/o.i  co2).  Temperature  (K) 


Fig.  9  illustrates  the  power-current  density  curve  of  the  cell 
for  different  ratios  of  carbon-dioxide  and  oxygen  (2.0, 1.0  and  0.5) 
fed  into  fuel  cell.  The  operating  temperature  and  pressure  were 
assumed  to  be  923  K  and  atmospheric,  respectively.  It  appears  that 
the  variation  of  this  ratio  has  no  significant  effect  on  the  cell  out¬ 
put  power.  However,  the  figure  demonstrates  a  slight  difference 
at  higher  ratios  which  indicates  that  the  excess  carbon  dioxide 
is  desired.  The  same  conditions  and  arguments  are  applicable  in 
describing  the  voltage-current  density  curve  (Fig.  10).  The  linear 
behaviour  of  the  polarization  curve  for  molten  carbonate  fuel  cells, 
as  is  proved  in  this  study,  is  normal.  The  reason  is  that  the  activation 
polarization  which  happens  at  low  current  densities  in  the  low- 
temperature  fuel  cell  (e.g.  proton  exchange  membrane  fuel  cell), 
disappears  in  a  high-temperature  fuel  cells.  This  observable  fact 
was  also  reported  in  some  other  studies  [13,17]. 

In  conclusion,  in  view  of  the  fact  that  the  current  density  is  a 
vector  quantity,  it  is  worth  to  mention  that  only  the  transverse 
component  is  valuable  and  contributes  to  the  output  power  of  the 
cell.  In  other  words,  the  lateral  component  only  diminishes  the 
output  power.  Precisely,  for  the  electrochemical  reaction  occurring 
in  the  triple  boundary  phase,  electrons  must  transfer  transversely 
through  the  electrode  to  be  collected.  Thus,  in  every  stage  of  molten 


Current  Density  (A  m'2) 


Fig.  10.  Effect  of  current  density  and  Pco2/Po2  ratio  on  cell  voltage  (operating  con¬ 
dition:  923  K,  fuel  gas:  0.7  H2/0.2  H2O/0.1  C02). 


Fig.  11.  Variations  of  cell  energy  and  exergy  efficiencies  at  different  operating  tem¬ 
peratures. 

carbonate  fuel  cell  design,  endeavour  must  be  implemented  to  min¬ 
imize  the  lateral  component. 

5.2.  Thermodynamic  results 

In  this  section,  a  number  of  parametric  studies  were  conducted 
in  details  with  the  intention  of  estimating  the  performance  of 
the  molten  carbonate  fuel  cell  with  different  operating  conditions. 
Using  equations  provided  earlier,  energy  and  exergy  efficiencies  of 
the  unit  cell  were  obtained.  In  addition,  the  entropy  generation  is 
investigated  in  some  cases.  In  general,  molar  fractions  of  hydrogen, 
water  vapour  and  carbon  dioxide  were  considered  to  be  0.72,  018 
and  0.1,  respectively.  Likewise,  oxygen  and  carbon  dioxide  molar 
fractions  were  chosen  to  be  0.33  and  0.67.  Anodic  and  cathodic  gas 
flow  rates  are  1.66  mol  h-1  and  2.04  mol  hr1,  respectively,  unless 
otherwise  stated. 

Variations  of  cell  energy  and  exergy  efficiencies  at  different 
operating  temperatures  (ranging  from  883  K  to  963  K)  are  pre¬ 
sented  in  Fig.  1 1 .  Results  were  obtained  based  on  operating  current 
density  of  4000 Am-2  and  atmospheric  pressure.  According  to 
Fig.  1 1 ,  energy  efficiency  of  the  unit  cell  varies  from  42.8%  to  50.5% 
while  the  exergy  efficiency  remains  in  the  range  of  26.8-36.3%. 
The  dissimilarity  of  the  two  aforementioned  efficiencies  is  caused 
by  internal  irreversibilities.  Perceptibly,  the  activation,  ohmic  and 
concentration  polarizations  are  the  major  source  of  irreversibili¬ 
ties  and  hence,  exergy  destroyed  inside  the  fuel  cell  assembly.  As 
a  result,  the  destroyed  exergy  diminishes  the  exergy  efficiency  of 
the  MCFC.  However,  it  is  apparent  that  both  energy  and  exergy  effi¬ 
ciencies  enhance  with  the  increase  of  temperature.  This  effect  can 
be  justified  with  the  fact  that  an  increase  in  operating  tempera¬ 
ture  reduces  the  irreversible  losses  (irreversibility)  of  the  fuel  cell, 
which  in  turn  results  in  augment  of  both  energy  and  exergy  efficien¬ 
cies.  Furthermore,  as  it  may  be  observed  from  Fig.  1 1,  even  though 
both  efficiencies  increase  with  temperature,  a  sharper  trend  takes 
place  in  lower  temperatures  and  they  tend  to  have  a  flatter  shape  at 
higher  temperatures.  In  fact,  the  unit  cell  total  losses  and  the  total 
input  energy  and  exergy  are  the  origin  of  this  alteration  in  the  effi¬ 
ciency  trends.  Clearly,  Eqs.  (29)-(31 )  illustrate  that  as  temperature 
increases,  the  activation,  ohmic  and  concentration  losses  (hence 
total  polarization)  drop.  The  higher  the  temperature  reaches,  the 
lower  the  irreversibilities  are.  In  contrast,  escalating  the  tempera- 
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Fig.  12.  Variations  of  cell  energy  and  exergy  efficiencies  at  different  current  densi¬ 
ties. 


sures. 


ture  results  in  a  constant  increase  of  energy  and  exergy  supplied  to 
the  unit  cell  while  the  net  power  output  of  the  cell  reaches  a  maxi¬ 
mum  and  then  declines.  As  a  consequence,  efficiencies  are  expected 
to  decrease.  There  always  seems  to  be  a  trade  off  between  the  effect 
of  ohmic  losses  and  supplied  energy  and  exergy. 

Fig.  1 2,  presents  the  effect  of  current  density  (ranging  from  1 000 
to  5000  Am-2)  on  energy  and  exergy  efficiencies.  For  this  case,  the 
operating  temperature  and  pressure  were  set  to  883  K  and  2  atm, 
respectively.  It  may  be  observed  that  both  energy  and  exergy  effi¬ 
ciencies  initially  raise  at  lower  current  densities  up  to  the  point  that 
they  attain  their  maximum  values  and  ultimately  decrease  with  the 
increase  of  current  density.  In  terms  of  energy  and  exergy  efficien¬ 
cies,  the  optimum  current  density  occurs  around  4000 Am-2.  At 
the  optimum  current  density,  58.7%  energy  efficiency  and  35.5% 
exergy  efficiency  were  achieved.  Since  the  operating  temperature 
is  considered  to  be  constant,  any  change  in  both  efficiencies  can  be 
attributed  to  the  net  power  output  of  the  MCFC  which  is  a  func¬ 
tion  of  cell  operating  voltage  and  current  density.  As  it  is  already 
discussed  in  Section  5.1,  the  voltage-current  density  curve  has  an 
increasing-decreasing  trend  which  concludes  a  similar  fashion  in 
both  efficiencies. 

Effects  of  operating  pressure  (1-5  atm)  on  energy  and  exergy 
efficiencies  of  the  studied  MCFC  are  illustrated  in  Fig.  13.  In  this 
parametric  study,  the  operating  temperature  and  current  density 
were  considered  to  be  923 K  and  4000 Am-2.  With  the  increase 
of  pressure,  both  the  energy  and  exergy  efficiencies  of  the  cell 
enhance.  This  is  realistic  since  as  pressure  increases,  an  extensive 
drop  in  irreversible  losses  occurs.  In  particular,  anode  and  cathode 
over-potentials  get  lower  values  in  higher  pressures.  To  be  more 
specific,  molar  concentration  of  hydrogen  in  anodic  triple  phase 
boundary,  and  oxygen  and  carbon  dioxide  in  cathodic  triple  phase 
boundary  increase  with  an  increase  in  operating  pressure.  As  a  con¬ 
sequence,  the  irreversible  losses  of  anode  and  cathode  decrease, 
which  in  turn  improve  the  performance  of  the  MCFC.  Neverthe¬ 
less,  both  efficiencies’  increase  is  sharper  at  the  operating  pressures 
below  3  atm  and  they  tend  to  have  a  smoother  change  above  this 
pressure. 

Fig.  14  describes  how  the  variation  of  anode/cathode  gas  flow 
ratio  can  affect  the  efficiencies  of  an  MCFC.  The  unit  cell  assumed 
to  operate  at  923  K,  atmospheric  pressure  and  4000  Am-2.  Results 
show  that  an  increase  in  this  stoichiometric  ratio  lessens  the  energy 
and  exergy  efficiencies  of  the  unit  cell.  These  changes  are  consider¬ 


Table  3 

Variations  of  entropy  generation  at  different  operating  temperatures  and  pressures. 


Temperature  (K) 

Entropy  generation 
(J  mol-1  K-1) 

Pressure  (atm) 

Entropy  generation 
(J  mol”1  K-1 ) 

883 

26.2 

1 

19.7 

903 

22.6 

2 

15.8 

923 

19.7 

3 

13.6 

943 

17.4 

4 

12.9 

963 

15.4 

5 

12.4 

able  for  energy  efficiency  rather  than  exergy  efficiency.  Moreover, 
it  may  be  observed  from  Fig.  14  that,  as  this  ratio  grows  over  the 
unity,  both  efficiencies  decrease  with  a  gentle  slope.  In  addition, 
both  efficiencies  have  their  peak  value  when  the  molar  flow  rate 
of  the  oxidant  entering  cathode  gas  channel  is  higher  than  the  fuel 
molar  flow  rate  which  gets  into  the  anode  gas  channel. 


Anode/Cathode  flow  rate  ratio 

Fig.  14.  Variations  of  cell  energy  and  exergy  efficiencies  at  different  anode/cathode 
flow  rate  ratios. 
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Table  3  summarizes  the  value  of  the  unit  cell  entropy  gen¬ 
eration  at  different  operating  temperatures  and  pressures.  As  it 
may  be  seen  from  this  table,  an  increase  in  operating  temperature 
decreases  the  entropy  generation  in  the  unit  cell.  This  effect  can 
be  acceptable  by  noting  the  fact  that,  the  entropy  generation  is  a 
function  of  cell  irreversibilities  in  terms  of  activation  and  ohmic 
polarization.  These  irreversibilities,  as  shown  in  Fig.  7,  are  consider¬ 
ably  reduced  with  the  increase  of  operating  temperature,  which  in 
turn  cause  a  drop  in  the  entropy  generation.  Likewise,  higher  oper¬ 
ating  pressures  result  in  lower  irreversibilities  and  hence  entropy 
generation. 

6.  Conclusions 

In  the  first  part  of  this  study,  a  mathematical  model  is  developed 
to  investigate  the  performance  characteristics  of  a  molten  carbon¬ 
ate  fuel  cell  at  very  high  current  densities,  by  employing  mass  and 
momentum  conservation,  electrochemical  reaction  mechanisms 
and  electric  charge  transfer.  A  finite  volume  based  commercial  soft¬ 
ware  package,  ANSYS  FLUENT  12.0.1,  is  used  to  solve  the  system 
of  PDE  equations.  In  addition,  the  model  was  applied  to  predict 
the  unit  cell  performance  characteristics  at  a  variety  of  operating 
pressures,  temperatures,  and  cathode  gas  stoichiometric  ratios.  The 
present  results  confirm  the  following: 

•  For  low  current  densities,  there  is  not  very  significant  difference 
in  output  voltage  and  power.  Flowever  for  greater  values  of  cur¬ 
rent  densities,  the  difference  is  not  negligible. 

•  The  fuel  cell  total  voltage  loss  drops  as  the  temperature  increases 
in  an  almost  linear  fashion  for  all  operating  current  densities. 

•  If  the  molten  carbonate  fuel  cell  is  to  operate  at  current  densi¬ 
ties  smaller  than  2500  Am-2,  there  is  no  point  to  pressurize  the 
system. 

•  Above  this  current  density,  the  output  power  is  enhanced  as  the 
operating  pressure  is  increased,  nevertheless. 

•  If  the  fuel  cell  operates  at  pressures  greater  than  atmospheric 
pressure,  the  unit  cell  cost  could  be  minimized. 

•  Various  partial  pressure  ratios  at  the  cathode  side,  demonstrated 
nearly  the  same  effect  on  the  performance  of  the  fuel  cell. 

In  the  second  part,  a  thermodynamic  model  is  utilized  to  examine 
energy  efficiency,  exergy  efficiency  and  entropy  generation  of  the 
MCFC.  Some  parametric  studies  are  performed  and  the  following 
results  are  obtained: 

•  By  changing  the  operating  temperature  from  883  K  to  963  K,  the 
energy  efficiency  of  the  unit  cell  varies  from  42.8%  to  50.5%  while 
the  exergy  efficiency  remains  in  the  range  of  26.8-36.3%. 

•  Both  efficiencies  initially  raise  at  lower  current  densities  up  to 
the  point  that  they  attain  their  maximum  values  and  ultimately 
decrease  with  the  increase  of  current  density. 

•  With  the  increase  of  pressure,  both  the  energy  and  exergy  effi¬ 
ciencies  of  the  cell  enhance. 

•  An  increase  in  this  anode/cathode  flow  ratio  lessens  the  energy 
and  exergy  efficiencies  of  the  unit  cell. 

•  Fligher  operating  pressure  and  temperature  decrease  the  unit  cell 
entropy  generation. 
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